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ABSTRACT 

1E2259+586 belongs to the gronp of anomalons X-ray pnlsars (AXP) which 
are thonght to be magnetars. Acenrate measnrements of the proper motion for 
these interesting objects have not been made. In this work we nse the data 
obtained by satellites ROSAT, Chandra and XMM-Aewton taken in a period 
of 12 years to measnre the change in position of 1E2259-I-586 with respect to 
the three held sonrees. Onr work yielded an npper limit on the proper motion 
of ~ 170 mas/yr which corresponds to a transverse speed of 2500 km s“^, for 
an assnmed distance of 3 kpc. We also nsed this npper limit to pnt a lower 
limit on the age of the CTB 109 of ~ 6 kyr, assnming they are associated. The 
consistency of right ascension and declination measnrements by different satellites 
snggest that we can average these measnrements to get the best to date position 
of 1E2259+586 . 

Subject headings: Anomalons X-ray Pnlsars: 1E2259-I-586 - AXPs: proper mo¬ 
tion - stars: magnetars 


1. Introduction 

Magnetars (AXPs & SGRs) being nentron stars, mnst be born dnring core-collapse 
snpernovae events. Hence a snpernova remnant is expected to be aronnd. Snre enongh 
1E2259-I-586 is associated with CTB109 and AXP 1845-0258 with SNR G29.6-1-0.1 (Woods 
& Thompson 2004). These associations assnme magnetar kick-velocities are similar to those 
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of radio pulsars. It would be very attractive to measure some magnetar proper motions 
directly. This work aims to do that. We have made use of the published and archival X-ray 
data to measure the proper motion of 1E2259-I-586. 

In order to eliminate pointing accuracy uncertainties of the satellites we measured the 
position of 1E2259-I-586 relative to other X-ray point-sources in the same held. In Table 1 
we list the archival X-ray observations where we could identify the right ascension (a) and 
declination (5) difference between 1E2259-I-586 and the three reference sources. 


2. Observations 

Since its hrst discovery Anomalous X-ray Pulsar (AXP) 1E2259-I-586 has been observed 
with various X-ray missions many times. These observations started as early as Einstein 
and continued all the way to state-of-the-art satellites like Chandra and XMM-iVewtou. For 
this work we choose the observations that satisfy the following criteria; 1) Data taken by 
detectors that have high spatial resolution (~ 10''), 2) All three reference sources (see Fig. 1) 
are in the held of view, 3) Exposure is long or the detector is sensitive enough so that all 
three reference sources are detected. In Table 1 we list the observations and the name of 
their respective observatories that we have used to get the positions of three selected sources 
together with the position of 1E2259-I-586. 

In order to get the positions of the selected sources in the held of view of the pulsar we 
extracted unbinned images of all hve observations. For data sets obtained with ROSAT PSPC 
and HRI, and Chandra ACIS we used the CIAO (version 3.2) tool WAVDETECT^ to detect 
all sources in the held of view. For XMM-Aewtou observation we used the SAS (version 
6.1.0) tool EBOXDETECT that performs a sliding box detection, using locally estimated 
background (see SAS data analysis threads web page^). Although using WAVDETECT on 
ROSAT HRI will not yield accurate Point Spread Function (PSF) sizes, this should not ehect 
our results because we are only interested in the central position of the sources. 


^http: / / asc.harvard.edu/ciao/threads/wavdetect / 

^http://xmm. vilspa.esa.es/external/xmm_sw_cal/sas_frame.shtml 
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3. Determining the Proper Motion of 1E2259+586 

In Table 1 we show the measured coordinates of all three reference sources together with 
the X-ray pulsar. Each of these measurements have been done for Eve different times. For 
each observation we calculated the distances of the reference sources to the pulsar in the di¬ 
rection of right ascension (a) and declination (h). The individual errors were calculated from 
the positional accuracies of the respective instruments. PSPC detector on board ROSAT 
has a spatial on-axis resolution of 25 arcsec FWHM at 0.93 keV. The other detector HRI, 
however, has relatively high spatial resolution of ~ 2". Later, more advanced instruments 
such as Chandra ACIS-S and XMM-Newton MOS have also good angular resolution. For 
an on-axis source ACIS-S point-spread function (PSF) is under sampled by the 0'.'492x0'.'492 
CCD pixels. On the other hand the EPIC-MOS camera has a spatial resolution of 5" FWHM, 
which is limited by the mirrors. Since all four sources have different spectra we chose these 
typical limiting resolutions as the error radii. These errors were then propagated with stan¬ 
dard error propagation to get the uncertainty on each angular separation. To each of these 
six data sets we made a error weighted least-squares fit with a straight line, y = ax + b (see 
Figure 2). 

Through the described method, for each data set we found the best fit value. The 
parameter a (slope) was varied around its best fit value (18 mas/yr) while b was kept constant 
at the value determined by the best fit. For each different value of a a value was calculated. 
This distribution is approximately a parabola. To calculate the combined distribution for 
all three data sets that give the proper motion for a (or 6) we then added the individual x^s 
and found the minima to get a best fit value (see Figure 3). To get the errors associated with 
the proper motions calculated we used the method described by Lampton et al. (1976). la 
error associated with the proper motion were found by adding 10.7^ to our minimum for 
total of 9 degrees of freedom and found the width of the curve at that point (Figure 3). We 
give the values of the best fit and it’s la error in Table 2. 


4. Absolute position of lE2259-|-586 

Since the data set on Table 1 represents also a summary of all the X-ray position mea¬ 
surements of 1E2259-I-586 , the consistency of the right ascension, declination measurements 
by different satellites suggest that we can average these measurements to get the best to 


^This is the value of the reduced chi-square ijC jv) corresponding to the probability of observing a value 
of chi-square larger than with v degrees of freedom. 
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date position of 1E2259+586 . The variance in the data set would also yield the error in 
this average position measurement. The results in (J2000) are, a = 23^01™08®21 ± 2'!2 and 
6 = +58°52'44''8±2'.'2. These coordinates, within the error-bars include candidates from pre¬ 
vious searches for the optical counterpart (Hulleman et al. (2001), Hulleman et al. (2000)). 
Optical identihcations, or upper limits are crucial for AXPs since it determines the mass and 
nature of the donor component, hence the energy emission mechanism. 


5. Discussion 

The uncertainties calculated being larger than the proper motion indicates that the 
available data is not sufficient to determine a dehnite proper motion. However, we can put 
an upper limit on it. The maximum angular speed 1E2259-I-586 can be calculated by the 
simple formula (yU^ -|- (see Table 2) which gives us an upper limit on the proper motion 

of ~ 170 mas/yr. This together with the adopted distance of 3 kpc (Kothes et al. 2002) 
gives us a maximum transverse velocity of 2500 km s“^. Although from the best fit values 
we can calculate a direction for the proper motion, the uncertainty being too high does not 
make the effort worth while. 

Associations, in general, between neutron stars and SNRs are usually judged on criteria 
such as agreement in age/distance, positional coincidence and evidence from proper motion. 
Distance estimates for AXPs have relatively big uncertainties, and there is no evidence 
that their characteristic ages are reliable age estimators. For 1E2259-I-586 we have two 
different age estimators; the age of the host remnant (~20 kyr, Rho & Petre (1997)) and 
also the spin-down age (~220 kyr). The age derived from the SNR assumes that they are 
physically associated. This method is a rough estimate for the age of the SNR and depends 
on the assumed spectral model for the X-ray emitting gas that surrounds the remnant. Our 
approach gives a lower limit for the age of 1E2259-I-586. Although we also assume that the 
SNR and the pulsar are associated, we make no assumptions on the spectral properties of 
CTB109 and its distance. The SNR is roughly spherical in shape and 1E2259-I-586 is close 
to its geometric center. We estimated the radius of this remnant as ~ 1000'' (see Figure 1). 
This radius together with the measured upper limit for the proper motion gives us a lower 
limit on the age of the pulsar of ~ 6 kyr. This limit is consistent with all other age estimates 
and is not very restrictive. However it is completely independent of all other methods used 
to determine the age. This limit is also very close to the age ~ 9 kyr determined by Sasaki 
et al. (2004) from several XMM-iVewtoii observations. 

The upper limit of 2500 km s“^ for the proper motion is consistent with that of those 
expected for radio pulsars. Arzoumanian et al. (2002) calculated the velocity distribution 
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of radio pulsars. They find that ~15% of pulsars have velocities greater than 1000 km s“^ 
and the distribution function asymptotically approaches zero beyond 700 km s“^ for a two 
component distribution. Although our upper limit is higher than all (but PSR B1552-31, 
(Manchester et al. 1978)) known pulsar velocities it is right on the upper limit of those that 
are expected for radio pulsars. 

Furthermore, birth properties of magnetars need not be the same as radio pulsars. For 
example, the magnetar model predicts that soft gamma repeaters may be born with velocities 
> 1000 km s“^ (e.g. see Thompson (2000)). We should also note that the velocity calculated 
here is a function of the distance, which is not well known and has an uncertainty > 50%. 

Another way of obtaining an estimate for the speed is to check the potential consequences 
high transversal velocities may cause. For instance, the speed of sound in the medium sur¬ 
rounding the pulsar is on the order of >5000 km s“^ assuming a gas temperature of 1 MK 
and mean particle density of 1 cm“^. If 1E2259-I-586 had a proper motion larger than this 
than we would expect a shock front to form, no such thing has been detected. Also, the pre¬ 
vious works that measured the absolute position of the pulsar for two different observations 
were not conclusive. The current best position (Patel et al. 2001) when compared to the 
previous best position (Hulleman et al. 2000) for 1E2259-I-586 yields no discrepancy within 
the uncertainties. Hence we conclude that the proper motion measurements that has been 
made in this work are the best so far. There are cases where optical sources are obscured by 
galactic dust absorbtion. Hard X-rays may be the only handle we have on the proper motion 
of some of these objects. However, relative proper motion measurements are usually done in 
optical and radio-astronomy since the requirements are: 1) there be enough sources in the 
field, 2) the accuracy of position measurements are high. The current observatories allows 
us to meet both of these criteria. Thus conhrming the coming of age in X-ray astronomy. 
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Fig. 1.— (Left)The region that is blown up here is shown with a box on the right hand side image. 
The sources around 1E2259+586 and their respective names adopted in this paper. The image is 
taken from Chandra observation, see Table 1. (Right) The circle approximately cocentered with 
SNR. The radius of the circle is 1000'', which is a rough estimate of the size of the remnant. 







Table 1: RA and DEC of selected sources around the pulsar, for different observations 


Mission 

ROSAT 

ROSAT 

Chandra 

XMM-JVewton 

XMM-JVewton 

Detector 

PSPC 

HRI 

ACIS 

MOS 

MOS 

Date 

1991-07-09 

1992-01-08 

2000-01-12 

2002-06-11 

2002-06-21 

Duration (ks) 

33.37 

22.41 

19.13 

52.53 

31.05 


Sources and their respective coordinates, (a,<5), for each observation 


1E2259-I-586 

345.2845 

-1-58.8783 

345.2828 

-1-58.8805 

345.2845 

+58.8792 

345.2842 

+58.8788 

345.2850 

+58.8788 

El 

345.2491 

-1-58.9628 

345.2498 

-1-58.9606 

345.2504 

+58.9589 

345.2492 

+58.9585 

345.2502 

+58.9586 

E2 

345.1378 

-1-58.8791 

345.1353 

-1-58.8815 

345.1386 

+58.8801 

345.1384 

+58.8797 

345.1397 

+58.8795 

E3 

345.1785 

-1-58.8416 

345.1784 

-1-58.8430 

345.1809 

+58.8426 

345.1804 

+58.8415 

345.1809 

+58.8414 
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Table 2: Statistics for the simultaneous fits done for Act and Ad for each selected source, see 
Figure 2. We used a linear function of the form y = ax + b. The errors given are la. 


xVdof 

a 

Aa 

Proper Motion {n) 


(arcsec/day) 

(arcsec/day) 

(mas/yr) 

ct 61.2/9 

-0.000054 

0.000285 

18 ± 104 

5 9.9/9 

-0.000033 

0.000285 

12 ± 104 
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Fig. 2.— The angular distance between 1E2259+586 and the respective X-ray source in the 
direction of a {Left) and 5 {Right) versus the date of the observation. The best linear fit to data 
is shown with a straight line in each case. We have taken out the datum for the ROSAT PSPC 
observation just for plotting purposes. It was, however, included in the fit. 
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Fig. 3.— psr degrees of freedom {u) versus proper motion in the direction of a {Left) and 
5 {Right). Solid curves are the distribution of the individual fits whereas the dotted line is 
the distribution of all three combined (average), is 3 for individual fits and 9 for the combined 
distribution. The dashed horizontal line is the marker for the la error level. See text for the 
explanation. 







